Strongly coupled large-angle stimulated Raman scattering of short laser pulse in 

plasma-filled capillary 



Serguei Kalmyko\Q 
Centre de Physique Theorique (UMR 7644 du CNRS), 
Ecole Polytechnique, 91128 Palaiseau cedex, France, and 
Max- Planck- Institut fur Quantenoptik, D-85748 Garching, Germany 

Patrick Mora 
Centre de Physique Theorique (UMR 7644 du CNRS), 
Ecole Polytechnique, 91128 Palaiseau cedex, France 
(Dated: February 2, 2008) 

Strongly coupled large-angle stimulated Raman scattering (LA SRS) of a short intense laser 
pulse develops in a plane plasma-filled capillary differently than in a plasma with open boundaries. 
Coupling the laser pulse to a capillary seeds the LA SRS in the forward direction (scattering angle 
smaller than 7r/2) and can thus produce a high instability level in the vicinity of the entrance plane. 
In addition, oblique mirror reflections off capillary walls partly suppress the lateral convection of 
scattered radiation and increase the growth rate of the SRS under arbitrary (not too small) angle. 
Hence, the saturated convective gain falls with an angle much slower than in an unbounded plasma 
and even for the near-forward SRS can be close to that of the direct backscatter. At a large distance, 
the LA SRS evolution in the interior of the capillary is dominated by quasi-one-dimensional leaky 
modes, whose damping is related to the leakage of scattered radiation through the walls. 

PACS numbers: 52.35 Mw, 52.38 Bv, 52.40 Fd 



I. INTRODUCTION 

The technique of chirped-pulse amplification pj made 
sub-picosecond laser pulses of high power (P > 10 12 W) 
available for generation of coherent x-rays , high har- 
monics of radiation (3L and laser wakefield acceleration 
(LWFA) of electrons 0, U |(j in rarefied plasmas [where 
luq ^S> <-<Jpe, is a laser frequency, uj pe = {Aire 2 riQ j 'm e ) x l' 1 
is an electron plasma frequency, no is a background 
electron density, m e and — |e| are the electron mass 
at rest and charge]. Full potential of these applica- 
tions can be realized with the laser-plasma interaction 
length increased beyond the Rayleigh diffraction length, 
Zr = 7t<7q/Ao, by means of external optical guiding jfj 
(here and later, A « 2ttc/ujo is a laser wavelength, and 
(To is a laser beam waist radius). In particular, a dielec- 
tric capillary, where the oblique mirror reflections sup- 
press the laser beam diffraction, can be used as a guiding 
tool 0, 0, 0] • Plasma can be created in a capillary by 
an optical field ionization of the filling gas [lfj, lllj, [EJ , 
or by a laser ablation of the walls ^3- Then, the large- 
angle stimulated Raman scattering (LA SRS) starts to 
challenge the transportation of a laser beam over a long 
distance pHT^. 

In the standard SRS process 0], the pump electro- 
magnetic wave (EMW) is scattered off spontaneous fluc- 
tuations of electron density, which, in turn, can be ampli- 
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fled by the ponderomotive beat wave of pump and scat- 
tered light. Appropriate phase matching of the waves 
results in a positive feedback loop with the onset of a 
spatio-temporal instability [T^. When the plasma ex- 
tent is much larger than a laser pulse length, and no 
reflections off plasma boundaries occur, both scattering 
electron plasma waves (EPW) and scattered EMW quit 
the region of amplification, and the convective gain sat- 
urates within a time interval of the order of pulse du- 
ration [H El HI E3 Given the pulse length, the 
maximum possible gain remains the same for all scat- 
tering angles, and whether it is achieved or not for a 
given an gle i s determined by the laser pulse aspect ratio 
only (l8l Il9t I20L l2l| . Convection of scattered radiation 
out of the laser waist may result in a strong pulse deple- 
tion [2^]. Even when the full depletion does not occur, 
the LA SRS can produce considerable pulse erosion [23| . 
suppression of the relativistic self-focusing |24| , heating 
and pre-acceleration of plasma electrons [25j , and seeding 
the forward SRS prl [23 • Thereby, knowing the details 
of the LA SRS evolution in various physical conditions is 
a matter of high importance for applications. 

Confining plasma between reflecting surfaces deeply 
modifies the SRS process. In the one-dimensional (ID) 
geometry, the Raman backscatter changes its nature from 
convective to absolute [2^]: reflections trap the unstable 
radiation modes inside plasma and give rise to the con- 
tinuous amplification. When the laser beam is confined 
between the mirror-reflecting partly transparent walls, 
and propagates collinearly to them, reflections reduce 
the sideward convection of scattered light. If the re- 
flective modes dominate in plasma, the LA SRS gain 
tends to that of the direct backscatter and thus re- 
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veals a dramatic increase in comparison with the open- 
boundary system The LA SRS in this geometry 
has been considered so far in the regime of weak cou- 
pling [l4|, when the scattering EPW is similar to the 
plasma natural mode [T?). and temporal growth rate is 
well below the electron plasma frequency. This regime 
requires fairly low amplitude of a laser pulse, i.e., clq <C 
y/ujp e /ujQ < 1 [ao = eEo j '(m e woc) is a normalized am- 
plitude of the laser electric field]. However, for the effi- 
cient LWFA the plasma density have to be re- 
duced in order to increase the gamma-factor of the laser 
wakefield, 7 g = LUo/ujp, and, hence, the electron energy 
gain. With l/j g < Oq < 1, the LA SRS becomes strongly 
coupled: its temporal growth rate exceeds aj pe , and the 
scattering EPW differs from the natural mode of plasma 
oscillations jlSL l2ll l2(il l29l |30( . In a plasma-filled cap- 
illary, the strongly coupled LA SRS acquires new spe- 
cific features: in a wide range of parameters relevant to 
the self-modulated LWFA y5j, our PIC simulations (us- 
ing the code WAKE |3l|) discovered a vast enhancement 
of the near-forward SRS in the immediate vicinity of the 
capillary entrance aperture. The unstable plasma modes 
were primarily transverse and therefore useless for the 
longitudinal electron acceleration. For the same range 
of parameters, this effect has never been significant in an 
open-boundary plasma. Independent fluid modelling [3^ | 
verified these observations. 

Making a step to understanding this phenomenon we 
propose a two-dimensional (2D) linear theory of strongly 
coupled SRS of a short laser pulse under a given angle a 
in a slab of rarefied plasma laterally confined between the 
mirror-reflecting partly transparent flat walls (flat capil- 
lary). Boundary conditions for the scattered radiation 
describe the oblique mirror reflections and the electro- 
magnetic (EM) seed at the entrance plane. We associate 
the latter with the signal formed of the high-order cap- 
illary eigenmodes produced by the laser beam coupling 
to the capillary (the coupling process is described else- 
where 0, E| an d is outlined in Appendix EJ . SRS in 
forward (a < 7r/2) and backward (a > 7r/2) direction 
proceed differently. Forward Raman amplification of the 
EM seed can be dominant within a finite distance from 
the entrance plane and be responsible for the instability 
enhancement observed in the modelling. On the other 
hand, the backward SRS is affected by the reflections 
only. As the LA SRS of a finite-length laser pulse pre- 
serves the convective nature (backward and, partly, side- 
ward convection of radiation is allowed), the gain satu- 
ration occurs within a finite distance from the entrance 
plane. Reflections give the unstable modes additional 
rise time in any transverse cross-section of a plasma, 
and, even for relatively small scattering angles (such as 
a = 7r/6 taken for numerical examples of this paper), 
the saturated convective gain can approach that of the 
backward SRS (BSRS). The field structure is then ap- 
proximated by a quasi-lD lossy mode whose damping is 
produced by the leakage of radiation through the walls. 

The paper is organized as follows. Section ITT1 presents 



a theoretical model of the strongly coupled LA SRS in 
a 2D slab geometry. The laser pulse entrance into a 
plasma and oblique mirror reflections of scattered light 
are expressed in terms of appropriate boundary-value 
conditions for the coupled- mode equations. General solu- 
tion of the boundary-value problem is presented (deriva- 
tion is given in Appendix [BJ| ■ Section II I II discusses the 
spatio-temporal evolution of instability. The case of fully 
transparent lateral boundaries is considered in subsec- 
tion IIII Al Enhancement of the LA SRS in the generic 
reflective case is considered, and the maximum gain fac- 
tors are evaluated in terms of appropriate asymptotic 
solutions in subsection IIII Bl Section Hvl summarizes the 
results. 



II. BASIC EQUATIONS AND SOLUTION OF 
BOUNDARY- VALUE PROBLEM 



In a 2D plasma slab confined between mirror-reflecting 
capillary walls, the high-frequency (hf) electric field is a 
superposition 



a(r,t) 



-iujQt 



+ C.C. 



(1) 



of that of the laser, ao, and of up- and down-going scat- 
tered EMWs, a s ± = eE s± /(m e ujQc) (|a s± | < |a | < 1); 
in Eq. and below, r = {x, z) is a radius- vector in a 
plane geometry. We assume that the waves Q have the 
linear polarization with the electric vectors parallel to 
the walls, and that the scattering occurs under an angle 
a in the plane orthogonal to the polarization vector. In 
the rarefied plasma, both ko = e z ko and k s ± obey the 
same dispersion relation loq = uj 2 e + crkh^ w c 2 k^ s± y, 
hence, |k a ±| = k s = k 0l k s±z = k sz = k cosa, ±k s ± x = 
k sx — fco since, and the amplitudes ao< s ±) vary slowly in 
time and space on the scales uJq 1 , k s ~ l , and kg^ 1 . Ions 
form a homogeneous positive background; this assump- 
tion holds for a laser pulse shorter than an ion plasma 
period, to -C 27rw pi 1 . The beat wave of incident and scat- 
tered radiation excites perturbations of electron density, 
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whose wave vectors obey the matching conditions k e ± = 
k -k s± ; hence, |k e± | = k e = 2k sin(a/2), k e±z = k ez = 
ko(l — cos a), k e ± x = ^fkosma. Wave vector diagram of 
the LA SRS is shown in Fig.^ The scattering EPW has a 
longitudinal component of phase velocity small compared 
to the speed of light, i.e., k ez > k p = u) pe /c. This restric- 
tion eliminates the forward Raman scatte ring |26ll33| and 
resonant modulational instability (RMI) |3J| . In the rar- 
efied plasma, the amplitudes N± vary slowly in space on 
the scales k e ~ l . kpZ 1 . 
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FIG. 1: Wave vector diagram of the LA SRS; "+"("-") stands 
for the up- (down-) going interaction. 



The amplitudes of up- and down-going scattered EMW 
and scattering EPW obey the coupled-mode equations 
derived from the equations of non-relativistic hydrody- 
namics of cold electron fluid in the hf field JU and 
Maxwell's equations for the scattered radiation, 



9 T, 9 , T, 9 

i I — + V Z — ±V X — | a s ± 
v ot, oz ox , 



+ kl\N± 



9iN±, (3a) 
g2d s ±, (3b) 



where V z — cosa/(l — cosa), V x — sina/(l — cosa), 
,9i = (a /2)(fc 2 /fc e J, and g 2 = a*(k e /2) 2 . The wave 
coupling parameter is G 3 = g\gi = (a /2) 2 k p k (strong 
coupling is the case for G ^> k p ). Equations © are 
expressed through the variables x, z, and £ = ct — z, that 
is, the temporal evolution of waves is traced in an x — y 
cross-section at a longitudinal position z. 

Figure |21 shows the interaction area. At z = 0, t = 0, 
the laser pulse enters a semi-infinite plasma-filled gap 
between flat mirror-reflecting walls (z > 0, < x < L x ) 
and propagates towards positive z. The pulse leading 
front, £ = 0, encounters the stationary level of electron 
density perturbations with a constant amplitude Nq, 



N±(x,z,Q) = N , 
dN±/d£(x,z,Q) = 0, 



(4a) 
(4b) 



fluctuations of radiation in fresh plasma being neglected, 



a s ±(x, z, 0) = 0. 



(5) 



At the capillary entrance plane, the transverse profile of 
radiation can have a significant content of the high-order 
capillary eigenmodes (coupling the incident laser beam to 
the capillary is discussed elsewhere 0, and is outlined 
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FIG. 2: Geometry of a laser pulse propagation in a laterally 
confined plasma. Laser pulse enters the plasma at z = and 
t = 0, and moves towards positive z. Boundary conditions 
are posed at the plasma boundary z = 0, the pulse leading 
front z = ct and the walls x = and x — L x . Rear edge of 
the pulse z = c(t — to) is a free boundary through which the 
waves quit the region of amplification. The boundary-value 
problem is solved in the area c(t — to) < z < ct, < x < L x . 



in Appendix . The resonant condition for the wave 
vectors selects the modes that can be amplified by the 
forward SRS (a < ir/2) thus leading to the formation of 
an EM seed signal. To facilitate the forthcoming analytic 
job, we give this signal in a simple parabolic form with an 
amplitude vanishing at the walls in order to provide the 
continuity of the solution in the interior of the capillary, 



a s ±{x, z = -0,£) = a S Q [l - (1 - 2x/L x ) 2 



(6) 



Inside the capillary, oblique mirror reflections couple up- 
and down-going EMW: each reflection converts an up- 
going wave into a down-going one and vice versa, 

a s +(0,z,O = r(a)a s -(0, z,£), (7a) 
a„-(L x ,z,£) = r{a)a s+ (L x ,z,£). (7b) 

The conditions Q set up a quasi-lD exponential behav- 
ior of waves at large z. The reflectivity coefficient is 
a known function of scattering angle, r{a) = |sina — 
l(5 w /5 pl f - cos 2 a] l ' 2 \/{sma + [{6 w /5 pl f - cos 2 a] 1 / 2 }, 
where S w and S p i are the refraction indexes of walls and 
plasma. Figure [3] shows r(a) for a glass capillary with 
S w rs 1.5 and S p i rs 1. 

The temporal increment of strongly coupled LA SRS 
exceeds the electron plasma frequency. Hence, we ne- 
glect k p N± in comparison with d 2 N±/d£ 2 in the left- 
hand side (LHS) of Eq. (|3b|) . With the allowance for 
not very tight capillary, the pump field envelope ao(x,£) 
represents a portion of laser radiation coupled to the 
capillary which experiences mostly paraxial propagation, 
k±/ko <C ^ujpe/too- Assuming that the pump field evo- 
lution at a given point (x, z) takes much longer than the 
SRS growth (zr/c 3> to), and in order to enable the an- 
alytic progress, we approximate ao(x, £) with a fixed flat 
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profile at any position z in a capillary of the width L x , below, H(y) is the Heaviside step-function. Solution of 
a o( x iO = o>qH[x)H(L x — x)H(£)H(cto — £), using the the boundary-value problem, 
effective pulse duration to and amplitude ciq. Here and 
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a s+ (R;r) = -a s0 JP(R; r) #a(; 1, 1/2; <0 - tJVo^ Y) — [1 - *ld(R; r, c*) - $ 2 zj(R; r, c,-)], (8a) 

1 N 3 

iV+(R;r) = - as052 ^(R;0(e-^/K) 2 oi ? 2(;2,3/2; J C) + ^^e c ^[l-$ 1D (R;r,c J )-3'2Z5(R;r,c,)], (8b) 

is then obtained via the 2D Laplace transform; here, R = (r, £), c 3 = iG 3 , C = (G 3 /4)(£ — z/V z ) 2 z/V z , 0-F2G &i,&2) 0) 
is the regularized generalized hypergeometric function |35j , 



- F a [ c ^Y^] \H{V z x-V x z) + Y,r n [H{V x z-V z x n . 1 )-H{V x z-V z , 



®2D = (l-r)£r n fi 



V x 



,Z)H(V x z-V z x n ), 



- , z \ f (V z x-V x z)(V z x-V x z-V z L x ) TT , T , Tr . 
^ = ^ HU -V z ) " few2? ^"^) 
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{V z L x /2f 



[H{V x z-V z x n - X )-H{V x z-V z x n )]y 
I 



(9b) 



(9c) 



where x n = x + nL Xl and the fundamental solution 
F s (jj,,v,£) (Ref. |22) is defined by Eq. (E7)l . The up- 
and down-going amplitudes are symmetric, a s _(x) = 
a s+ (L a; — x), iV_ (.t) = N + {L X — x), so we consider below 
the evolution up-going waves only. Our solution pos- 
sesses the same generic structure as the weakly coupled 
reflective solution discussed in detail in Ref.[l4|. However, 
contrary to the case of semi- infinite laser pulse of Ref.lbH 
the growth time of unstable waves is now limited by the 
pulse duration to, and, as shown in the subsection 1111 Bl 
the gain at a given point (x, z) remains finite and is given 
by either Eq. (JT3JI or JT§J|. 



III. SPATIO-TEMPORAL EVOLUTION OF 
UNSTABLE WAVES 



A. Evolution of instability in open-boundary 

system 



When the lateral boundaries are fully transparent, r = 
0, the EM seed at the entrance plane vanishes (a S Q = 0), 
and the instability grows from the electron density noise 
in a fresh plasma ahead of the pulse. The functions @ 



then read 

<P 1D (R;r, Cj ) = F s ( Cj ,z/V z ,Z)H(V z x-V x z), 
$ 223 (R;r, Ci ) = F s ( Cj ,x/V x ,£)H(V x z -V z x). 

Plasma is divided by the characteristics £ = z/V z , x = 
z(V x /V z ), £ = x/V x , into ranges of dependence, where 
the solution is prescribed by the boundary conditions for 
radiation posed at the pulse leading edge £ = (range I), 



£ < x/V x 
£ < z/V z 
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0, $ 



2D 



at the wall x — (range II), 
x < z(V x /V z ) \ 

£ > x/V x = °' $ 2D # 

or at the entrance plane z = — (range III), 



(10) 



(11) 



(12) 



The principal feature that makes the LA SRS of a short 
pulse |2lJ different from the case of semi-infinite laser 
beam |l9j is the gain saturation within a finite distance 
from the entrance plane. At some point, z sa t < +00, the 
scattered radiation arriving from the plasma boundary 
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FIG. 3: Coefficient of reflection versus scattering angle for a 
glass wall with an index of refraction 5^, = 1.5 (the refraction 
index of plasma is taken equal to unity). 

z = — drops behind the laser pulse, and in all the points 
z > z sa t neither part of the pulse belongs to the range III. 
Then, the evolution of waves and the gain do not alter 
with z. The gain saturates differently for the forward 
(a < ir/2) and backward (a > tt/2) scattering. 

When a < ir/2, and the distance from the entrance 
plane is not too large, i.e., z < min {V z cto, L X (V Z /V X )}, 
all the three areas (llL)fl - <|12[) are available in the pulse 
body (that is, within a rectangle < £, < cto, < x < 
L x ). Given the point (x,z), waves fall initially within a 
range I, where they grow in time exponentially with an 
angle-independent increment 

70 = (V3c/2)G « (V3/2) (a a /2) 2 ^ e . (13) 

Note, that k = 7o/c is the known "spatial" incre- 
ment of the strongly coupled BSRS in the co-moving 
frame [E H |H M, HI 113 The evolution of waves 
is strictly ID in space on this stage. Later, informa- 
tion from the boundaries x = and z — — reaches the 
point (x,z), and the spatial dependence becomes either 
2D for £ > x/V x , x < z(V x /V z ) (range II) or remains 
ID for £ > z/Vg, x > z(V x /V z ) (range III). The waves 
are not exponentially growing at this time. In the range 
III, the entrance effect dominates: vanishing the scat- 
tered EMW at the boundary z = — determines the 
behavior of ID amplitudes. Deeply enough in plasma, 
z > min {V z cto, L X (V Z /V X )}, the entrance effect vanishes 
as the pulse terminates sooner than the scattered EMW 
from the entrance plane can reach the observer at a given 
z. The pulse body is then divided between the ranges I 
and II, and the evolution of LA SRS is the same through 
the rest of the plasma. 

For a > 7r/2, the boundary- value condition posed for 
radiation at z = — can only produce the scattered EMW 
convecting outwards (the EMW characteristic £ = z/V z 
recasts in the lab frame variables as z = —ct \ cosa|, and 
corresponds to the wave propagating towards negative z). 
Thus, the entrance effect does not change the solution 
at a positive z, and the boundary- value condition © 



becomes excessive (this is also valid in the reflective case). 
The spatio-temporal evolution of the LA SRS remains the 
same at any z > 0. 

Given the pulse duration, the maximum of the sat- 
urated gain, a s+ ,N + ~ e 70 * , does not depend on the 
scattering angle, and whether it is achieved or not is de- 
termined solely by the pulse aspect ratio [2lJ . If the pulse 
is wide, or the scattering angle is sufficiently large, a > 
ctQ = 2 arctan (cto/L x ), the maximum gain is achieved at 
the pulse rear edge £ = ct for ctocot(a/2) < x < L x . 
Hence, the angular spectrum of scattered light is pre- 
scribed by the pulse aspect ratio rather than the angular 
dependence of the increment. For L x cto, scattering 
within a broad range of angles 2cio / L x < a < tt proceeds 
with the maximum gain. Otherwise, for L x <C cto, the 
highest gain corresponds to the near-backward scattering 
only, 7r — L x /(cto) < a < it (Ref. \2w . To estimate the 
pulse energy depletion due to the LA SRS, it is sufficient 
to neglect the radiation scattered under angles smaller 
than a fRef.l36h. 

B. LA SRS evolution in the reflective case 

In a capillary, the oblique mirror reflections of scat- 
tered light off the walls [the boundary condition 0] con- 
tribute to the LA SRS evolution over the whole laser 
path in plasma, but become actually dominating later, 
when the entrance effect vanishes (z > V z cto). The re- 
flections establish a long-distance asymptotic state of the 
LA SRS — amplification of a quasi-lD radiation and 
plasma modes with a temporal increment close to that 
of the direct backscatter l|13l) . Besides, the forward SRS 
is seeded by the EM signal a s+ (x,— 0,£) at the capil- 
lary entrance plane [the boundary condition ijojl]. If the 
signal amplitude a s o exceeds the amplitude of the elec- 
tron density noise No, the unstable waves will achieve a 
large amplitude (by virtue of the high seed level) at the 
transient stage of laser propagation, < z < V z cto- In 
this case, contrary to the SRS in the unbounded plasmas 
described in the preceding subsection, the near-forward 
SRS exhibits much higher level of amplification than the 
backward SRS (a > ir/2), and its contribution to the 
dynamics of electron density perturbations can be dom- 
inating. Nonlinearities of plasma response that can then 
appear are worth investigating and will be addressed in 
future publications. 

The asymptotic 0-^2 (; bi, b 2 ; 9) ~ 

(2 7 rx/3l_ 1 e( 1 -^-^/ 3 e 3 ^ + O(0-^ 3 ) at |0| -> 00 
(Ref. 1351) helps to evaluate the level of the EM seed 
amplification within the interval < z < V z cto, 

T e 3V3C 1/3 /2 
\a s+ (z < V z ct )\ - a s0^= — ^ — , (I 4 ) 

where C » 1. Then, \N+\ ~ {V z / gi )(C 1/3 /z)\a s+ \. At 
the pulse trailing edge, £0 — cto : the argument of the 
asymptotic reaches the maximum £ max = (Gcto/3) 3 at 
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Zmax = V z cto/3. The EPW amplitude at this point is 

-™+(?0,^max) ~ a S 7= ^-7- , (15) 

and |a s +(^o,2max)| ~ (gi/G)\N+(£ , z max )\. Equa- 
tion 1)151) shows that the maximum gain on the tran- 
sient stage is determined by the ID temporal incre- 
ment 1)1 3)1 independent on a (the angular dependence is 
retained in the pre-exponential factors only). For the 
minimally allowed scattering angle, a m i n = yj2us pe /(jjQ 
(hence, V z « wo/wp e ), the point of the maximal gain is 
Zmax(o:min) ~ (^o A^pe) (eta /3) 3> cto- Parameters of the 
following numerical example (Fig. give 2 m ax(«min) ~ 
3 mm, which is shorter than a typical ca pillary length 
used in experiments (> 1 cm) 0, S EH Cul ll^, E 01 ■ 



For z > 



the entrance effect becomes less pro- 



nounced and finally vanishes at the point z = V z cto — 
3z max , where the EM signal arriving from z = — drops 
behind the laser pulse, and the instability growth sat- 
urates. Evolution of both forward and backward scat- 
tering is then determined by the lateral reflections only 
and, given the scattering angle, remains the same in any 
x — y cross-section for z > V z cto- We show in Appendix ICl 
that at z > V z cto, and 70^0 > 1, a cumbersome exact so- 
lution (|8a|l tends asymptotically to a quasi-lD damped 
mode with a temporal growth rate close to 1)13)1 that gives 
the asymptote for the density perturbation amplitude, 

N+(x,£) ~ {No/3)[(l-r)/lnr]e s ^- lnr ^ L "\ (16) 

Here, s = (70 - Aj)/c, Aj/c » -V x \nr/{3L X ). Equa- 
tion (|16fl is valid under the "low leakage" condition, 
r > exp(—3GL x /V x ), which indicates that the scattered 
light is mostly trapped inside a plasma slab: the energy 
leakage through the wall at one reflection (that produces 
an effective decrement A7 <C 70) is less than the energy 
gain on the way between the walls. The lateral growth of 
the asymptote (|16|) is much slower than the growth with 
£. Asymptotic solution 1)16)1 displays the basic result of 
the reflective theory of the LA SRS: at large distances 
from the entrance plane and large coefficients of amplifi- 
cation, the amplitudes of unstable waves tend to quasi- 
lD leaky modes exponentially growing in time with the 
increment tending to that of the BSRS 1)13)1 . 

Comparison of Eqs. 1)15)1 and 1)16)1 shows that the max- 
imum amplitude of the scattering EPW on the tran- 
sient stage, z < V z cto, differs from the final asymp- 
totic level of density perturbations roughly by a fac- 
tor of (G/5 , i)(7oio)~ 1 ^ 2 (aso/A^o)- When this factor is 
larger than unity, and \N + (z > V z cto,x,cto)\ ~ 1, the 
plasma response can become nonlinear in the vicinity of 
z = V z ct /3. This could be avoided by keeping the ra- 
tio of the seed amplitudes a s o/N below (ffi/G)(7oto) 1 ^ 2 • 
The amplitude of the electron density noise is difficult to 
control in experiment; however, as shown in AppcndixlAl 
the content of the high-order eigenmodes in the laser radi- 
ation coupled to the capillary (and, hence, the amplitude 
a S Q of the EM seed signal), can be effectively reduced by 
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FIG. 4: Spatio-temporal evolution of the up-going EPW in 
the field of transversely limited laser pulse of finite duration; 
the pulse aspect ratio is cto/L x — 0.5. The scattering angle 
is q = 7r/6. Temporal evolution of the amplitudes is traced 
at the longitudinal positions (a) z = VzCto/3 « 2.2cto, (c), 
(d) z = V x cto « 6.5cto- The left column shows the SRS evo- 
lution inside the glass capillary (r = 0.42) with the EM seed 
amplitude at the entrance plane a s o = 0.4 x 10~ 5 . Plot (b) 
shows the saturated solution for an open-boundary plasma, 
r = o s o = 0, for z > (V z /V x )cto w 3.4cto- In the plot (d), the 
long-scale asymptote 1161 of the reflective problem (dashed 
line) is compared with the exact solution (solid line) and with 
non-reflective (dotted line) and BSRS (dash-dotted line) so- 
lutions at the point z — V z cto, x = L x /4. 



increasing the capillary radius versus the radius of the 
incident laser beam. 

Figure 21 shows the spatio-temporal evolution of an 
up-going EPW [Eq. l|8b"|)] for the SRS under the angle 
a = it/ 6. The laser and plasma parameters are ao = 0.7, 
Ao = 0.5 /Jin, cto = 0-5L x — 6fc~ 1 , u! pe /u>o — 0.007, which 
2.2 x 10 17 cm" 3 , the pulse duration t w 230 fs, 



give n 

and the maximum increment 70 ~ 2.25uj pe . The level of 
EM seed is chosen a s o ~ 0.58 x 10~ 5 a (according to Ap- 
pendix^ it corresponds to a capillary by a factor of two 
wider than in the case of perfect matching; by the defi- 
nition, in axi-symmetric geometry, the perfect matching 
condition provides coupling 98% of energy of an incident 
Gaussian laser pulse to the fundamental eigenmode EHn 
of a capillary y|). The level of the plasma noise evalu- 
ated in Appendix El is N ~ 1.5 x 10 -6 . The plots (a) 
and (c) correspond to the plasma confined in a glass cap- 
illary with the reflection coefficient r = 0.42 (see Fig. |3J, 
and (b) — to the unbound plasma (r — and a s o = 0); 
plot (d) shows the long-term asymptotic behavior of the 
reflective solution. The plasma cross-sections are set at 
(a) z = V z ct /3 w 2.2cf , (c), (d) z = V z ct w 6.5ct - 
Given the calculation parameters, the non-reflective so- 
lution saturates at z = L X (V Z /V X ) « 3.4c£q and is exactly 
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FIG. 5: Ranges of dependence for the SRS under the angle 
a = 7r/6 in the cross-sections at (a) z = V z cto/3 ~ 2.2cio and 
(b) z = V x ct w 6.5ct . The thick solid line x = z(V x /Vz) 
divides the principal ranges prescribed by the non-reflective 
theory. The dashed lines are the characteristics of reflected 
waves, x + nL x = z(V x /V z ), £ = z/V g , and £ = (a; + nL x )/K,. 
One reflection contributes to the scattering process in the 
range II in the case (a), and two reflections — in the case (b). 
Raman amplification of the EM signal given at the boundary 
z = —0 contributes to the solution in the sub-range lib and 
range III. 



the same in any plasma cross-section x — y beyond that 
point. Figure Efb) shows this solution. The ranges of 
influence of boundary conditions are shown in Fig. [SJ In 
the range I [see Eq. <jTU|) ]. the electron density noise from 
£ = is amplified, and the waves do not experience re- 
flections. In the range Ila the instability is seeded by the 
free-plasma noise, and yet is enhanced by the reflections; 
the Raman amplified signal arriving from the entrance 
plane is added to these waves in the range lib. The EM 
seed from the entrance plane z = — is amplified by the 
forward SRS in the range III. 

In a capillary, the forward SRS is considerably en- 
hanced at z < V z cto (roughly by a factor of 70 in 
amplitude) versus the case of unbound plasma [com- 
pare Figs. Ufa) and (b)]. Despite the EM seed am- 
plified in the range lib is non-exponcntially growing, 
the high ratio of seed amplitudes, a s0 /N ss 2.7 » 
(gi/G)(7oio) 1 ^ 2 ~ 0.026, makes the entrance effect rather 
pronounced [plot Ufa)]. For N w 1.5 x 10~ 6 , and T ~ 1, 
Eq. JTSJl gives log 10 \N+(£ , z max )/N \ w 7.15, which 
agrees with Fig. Ufa) . Hence, parameters of the numer- 
ical example lay at the border of validity of the linear 
approach, and increase in a s o will result in the nonlincar- 
ity of the plasma response on the transient stage {e.g., 
perfect matching gives ao s ~ 0.012ao, hence, according 
to Eq. (|n>J), \N + (£, ,z max )\ ~ 10 6 ; this burst of the for- 
ward SRS has been the regular feature in our numerical 
experiments and in the fluid simulations H^]}. On the 



other hand, reduction a s o to the level 7 x 10~ 8 ao (which 
in the axi-symmetric case would correspond to a capil- 
lary tube by a factor 2.5 wider than in the case of perfect 
matching, see Appendix^), will make the forward Ra- 
man amplification of the EM seed almost negligible and 
thus tolerable on the transient stage. 

Plot^c) shows a quasi- ID saturated solution [compare 
with Fig.^Jb)] thus shaped by the contribution from two 
reflections [according to Fig. E^b)] , which, in full agree- 
ment with the long-scale asymptote (|16H . demonstrates 
the growth rate close to that of BSRS. Difference between 
Figs. Ufa) and (c) shows that, under the parameters of 
our example, the forward scattering (a < 7t/2) is charac- 
terized at z < V z cto by much higher gain than the SRS 
in the backward direction (a > 7r/2). This situation is 
completely reverse of the SRS in an unbounded plasma, 
where the gain can only fall as an angle drops [Ig. So, 
the higher the EM seed level produced by the laser beam 
coupling (that is, the tighter the capillary in a numerical 
or real-scale experiment), the more important becomes 
the forward SRS. In such case, a high amplification level 
of waves may be observed within quite a long distance 
in plasma, z < V z (a min )cto « (u)q /u> pe )cto [see also the 
discussion following Eq. (|15[) ]. This effect is adverse for 
such applications as the self-modulated LWFA in capil- 
laries 15] . The way of reducing the excessive forward 
SRS enhancement can be found in using a wider capil- 
lary (the laser focal spot fixed) than the perfect matching 
requires. 

Figure Old) traces the temporal evolution of the up- 
going EPW at x = L x /A (near the capillary wall). 
The asymptote ijTBj) perfectly approximates (and, for 
applications, can be used instead of) the exact reflec- 
tive solution (|8 b|) ■ The dash-dotted and dotted lines in 
Fig. gfd) correspond to the BSRS solution |AT + (£)| w 
(No/ 3) exp (7o£/c) and the exact non-reflective solution, 
respectively, providing the upper and lower limits of the 
convective gain variation. In this example, only in the 
vicinity of the border x — the coefficients of amplifi- 
cation in the reflective and non-reflective cases are con- 
siderably different [compare Figs. Hfo) and (c)]. Contri- 
bution from the reflections increases the wave amplitude 
at x = L x /4 by roughly an order of magnitude (compare 
solid and dotted lines at k p t; w 6). For the parameters 
chosen, the scattering EPW remains linear in the con- 
vective saturated regime (z > V z cto): the plasma noise 
level, 7V w 1.5 x 10~ 6 , substituted into Eq. (|Pd|) . gives 
\N + \ < 0.2 throughout the whole time interval < t < t . 

In the limit r — > 1 the total suppression of the lat- 
eral convection occurs. The up- and down-going ampli- 
tudes ijSJ) become purely one-dimensional for z > V z cto 
and completely identical. These amplitudes grow in time 
exponentially with the BSRS increment l|13|) . 
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IV. CONCLUSION 

We have proposed a 2D non-stationary linear theory 
of strongly coupled LA SRS of a short laser pulse in a 
flat plasma slab confined between mirror-reflecting walls 
(flat capillary). In a capillary the lateral convection of 
scattered light is partly suppressed by the oblique reflec- 
tions, and the instability experiences an enhancement. 
Additional enhancement of the SRS in forward direction 
(a > 7r/2) is produced by the amplification of the elec- 
tromagnetic seed signal that is formed of the high-order 
capillary eigenmodes at the entrance plane (formation of 
the signal is a consequence of the laser beam coupling to 
the capillary). The convective nature of LA SRS does not 
change. The asymptotic behavior of the waves demon- 
strates the transition from the set of 2D modes to the 
dominant quasi-lD damped mode. Even for near-forward 
scattering the convective gain of the dominant quasi- ID 
mode may be close to the BSRS gain. 
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APPENDIX A: SEED SOURCES FOR LA SRS 

The LA SRS under arbitrary angle in a strongly rar- 
efied plasmas (aj pe <C wo) is seeded by spontaneous elec- 
tron density fluctuations ahead of the pulse. A root- 
mean-square (rms) amplitude of these fluctuations is rep- 
resented in the equations by the quantity No, which 
gives the amount of seed corresponding to the element 
of solid angle dSlk,. in the direction of the wave vec- 
tor k e of scattering EPW, and can be expressed as 
N 2 {k e ) = dn kc Jn 3 e D k 2 dk sa n 3 e D (k e )k 2 e Ak e dn kc , where 
n 3D (k e ) is a spectral densi ty o f electron fluctuations in- 
tegrated over frequencies (331 • The integral is taken 
over the area of maximal spectral density of scattering 
EPW [k p=> k e , the amplification bandwidth cAfc e ps 
4v7o/^o is estimated at 7oio 3> 1 with taking account 
of the gain narrowing (Ref. Isfih], Using the ratio of 
the phase volumes \dk e /dk s \ = 2 sin 2 (a/2) (Ref. l38lh 
we express the seed amplitude through the element of 
solid angle c?f2k a in the direction of detector, Nq ps 
8nl D (k e )k%[(7o/c)/{ch)] 1/2 sm 2 (a/2)dri ks . Our theoret- 
ical formalism based on the assumption of quasi- plane in- 
teracting waves requires small variation of the scattered 
wave amplitude across the direction k s in the transversely 
limited area, < x < L x , which gives an estimate of 
the angular spread Aq ps cos a/(L x ko). The element 
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FIG. 6: Overlap integral between the radial profiles of am- 
plitude of the incident laser (Gaussian) and of the capillary 
hybrid eigenmodes EHi„. 

of solid angle then estimated as A£!k s ~ 2tt sin aAa — 
7r sin 2a/(L x k ) gives 

/ 8^ l + {k e r De ) 2 { l0 t y/ 2 sin 2 (a/2)sin2a 
1 0| ~Y n X 3 Q 2+(k e r De ) 2 Mo) 3 / 2 k L x 

( A1 ) 

where the spectral density of low-frequency electron fluc- 
tuations n 3D (k e ) is evaluated using formula (11.2.6.6) 
of Ref. [33- Parameters of Fig. 0] and k e ro e <§C 1 give 
|JV | fa 1.5 x 1CT 6 . 

For the strongly coupled SRS in the forward direction 
(a < 7r/2), coupling the laser beam to a capillary creates 
additional source of instability. The radial profile of an 
incident radiation with the wings cut off by the edges of 
the entrance aperture is approximated with an expansion 
through an infinite number of radial eigenmodes (having 
the same frequency wo) BUS- The high-order eigenmodes 
(characterized by the frequency u> 8 — ujq and high trans- 
verse wavenumbers, k n ± ~ k sx , where the integer n is 
the mode order) form the seed signal that is further am- 
plified in plasma by the SRS [in the model form, the 
transverse profile of this signal is given by Eq. ©]. It 
should be emphasized that these modes provide no seed 
for the weakly coupled SRS (including near-forward SRS 
and RMI [34[ that correspond to small scattering angles, 
a <C ^/ujpe/ujo), as this process requires the frequency 
matching uj s ps luq — w pe between the seed and the pump. 

Exact functional form of the capillary eigenmodes de- 
pends on the geometry chosen. Despite the LA SRS in 
a flat capillary is considered in the paper, we suppose 
that an estimate of the EM seed level will be more use- 
ful for applications if inferred from the axi-symmetric 
theory of the laser beam propagation in a dielectric 
tube 0, II] • The theory represents an electric field pro- 
file at the entrance aperture of the radius tq as an in- 
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finite sum a(r) — a J2^ 1 C n jQ(k n ±r) fRef. fTTf> . where 
C n = 2[r J 1 (u n )]~ 2 J r ° a(r)J (u n r/r )rdr is the overlap 
integral of the hybrid capillary eigenmode EHi„ with the 
incident laser profile a(r) = exp(— r 2 jofy (Fig. 0). Here, 
u n is the nth zero of the zero-order Bessel function of 
the first kind, Jo(u n ) = 0; for k r ^> 1, and n ^> 1, 
k n ±ro ~ u n m (n + 1/2)tt. Fig. [|)]shows that about 98% 
of laser energy is coupled to the fundamental mode for 
(Jo = 0.645ro (the perfect matching condition); however, 
the overlap integral decays very slowly as n grows {the 
analytic fit, C n ~ F n = exp[— (n + 50) 3 ], is almost ex- 
act for n > 10}. When the ratio a /r drops, the larger 
number of lower-order eigenmodes is effectively excited 
(up to 5 for a = 0.258r ), but contribution from the 
higher-order modes into the radiation profile at z — — 
drops sharply (e.g., analytic fit C n ~ 3 x 10 _6 F„, n > 10, 
holds for (To = 0.258ro). Therefore, choosing wider capil- 
lary is the way of reducing the effect of laser coupling on 
the forward strongly coupled SRS. 

The numerical example of subsection IIII Bl shows how 
the SRS under the angle a — ir/6 develops in the cap- 
illary with an EM seed amplitude characteristic of the 
axi-symmetric capillary tube by a factor of two wider 
than in the case of perfect matching. For the parameters 
of Fig.0J with the value L x /2 SdOk^ 1 assigned to r , 
equalizing the effective radial wave number k n ± ~ nir/ro 
to the resonant wave number k s ± = fco sin a gives the 
resonant mode order n* sa (fco?"o/7r) sin a 136 cor- 
responding to C136 ~ 0.58 x 10~ 5 under the condi- 
tion (T = 0.3225r (analytic fit C„ w 0.7 x 10~ 3 F„ 
was used; see also Fig. [SJ). One can expect that sev- 



eral modes with n n* can contribute to the effec- 
tive amplitude a s0 of the seed signal. The difference 
in the mode numbers An = n* — n causes the angular 
spread Aa around the given scattering angle; this spread 
should not exceed the admissible value established above, 
Aa < cosa/ (L x k ) ~ cos a/(2r ko). For |Aa| <C a 
and |An| <C n* , one has An ~ (fco r o/ 7r ) cosaAa < 
cos 2 a/(27r) < 1, so that only one capillary mode with 
n = n* can contribute to the scattering under the given 
angle. Therefore, the seed amplitude used for the nu- 
merical demonstration of Fig. QJa) is evaluated as a s o ~ 
C 136 a r> 0.4 x 10- 5 . 



APPENDIX B: DERIVATION OF THE EXACT 
REFLECTIVE SOLUTION 



Omitted k 2 in the LHS of Eq. (|3b|) , the Laplace trans- 
form (LT) of Eqs. © with respect to £ (LT variable s) 
and with respect to z (LT variable p) gives the set of 
ODE for the Laplace images a s ±(x,p, s; r), 



d_ 

dx 



T- )a s ± = ±K(x), 



(Bl) 



where il = (V Z /V X )(T S - p), T s = (iG 3 /s 2 - s)/V z , 
K(x) = Kj/isp) + (K 2 /s)[l - (2x/L x - l) 2 ], K x = 
giN /(iV x ), K 2 = (V z /V x )a s0 . The boundary condi- 
tions are a s+ (0,p, s) = ra s _(0,p, s) and a s -(L x ,p, s) — 
ra s+ (L x ,p, s). Equations IjBip admit the solution 



J 



a s+ {x,p,s;r) = - [ — + K 2 
P 



1 - 



(l-r)e &1 
1 — re nLi 



{L x /2) 2 \-re nL A \ 2 




(B2) 



where x n = x + nL x . Lateral symmetry gives a s -(x) = a s+ (L x — x). The Laplace transform inversion of Eq. I|B2 
with respect to p reads 



a s+ {x,z,s;r) = — — - 
si s 



K. — -K2 
si .« s 



~ e lsZ (x-z)(x-z-L x ) 



(L x /2) 2 



H(x-z)+Ki(l-r)^Tr n 



,r a (v z /v x )xr, 



n=0 



sT s 



-H(z-Xn) 



K 2 — {Xn ~\^ 5) h ~ Z) \ [H(z - x n . x ) H(z x n )] 



sT 



{L x /2f 



(B3) 



where, in accordance with Ref. Il4i the expansion [1 — 
r exp(riL a; )] _1 = ^^L r" exp(nf2La;) is used, and z = 

(V x /V z )z, K h2 = (V x /V z )K h2 - Equation includes 
Laplace images of the three types: l/(sr s ), e r ^/(sr s ), 



and e TsV /s. Their inversions read: 

CJ 1 {e r °y/s} = V^H(£ - z/V z ) F 2 (- 1, 1/2; »0, (B4) 

expressed through the regularized generalized hyper- 
geometric function [3^ of variable i£ = i(G 3 /4)(^ — 
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z/V z fz/V z ; 




where c| = iG 3 , so that c\ = -iG, c 2 = (i + V3)G/2, 
c 3 = (i - V3)G/2; and 




which is expressed through the fundamental solutions |2 



F s ( M)U) 0=e-^^-^- 7 (2n+l,^-,;))^-«), 

n=0 

(B7) 

where 7(771, r) = J e T t' dr 1 is the incomplete 
gamma-function of the order m of a complex variable 
t (Ref. 35). Combining expressions (|B5fl . (|B6|) . and l|B4|) 
in Eq. I|B3(I gives the envelope (|8a|l of the up-going EMW. 
The amplitude l|8b(l of the scattering EPW is derived in 
the similar fashion. 



APPENDIX C: ASYMPTOTIC REFLECTIVE 
SOLUTION 

The asymptotic can be found by applying the inversion 
formula, a s+ (x,z,s) = (27rz)~ 1 f°_f™ e px a s+ (x, p, s) dp, 
to the expression i|B2|l . The asymptotic behavior at 
z — > 00 is determined by the singularities of the integrand 
at p = 0, p — r s , and p n = T s - v + (V x /V z ){2imi/ L x ) 
with n integer and v = — (V x /V z )(hir/L x ). Expand- 
ing A s+ (x,p, s) in the vicinity of the specific points 
shows that all the singularities at p — T s and p = p n 
give the contribution ai(x,z,s) cx exp(r s z). Therefore, 
£j 1 {ai} cx H(£ — z/V z ), and, for a pulse of finite dura- 
tion, £ < cto, no contribution comes to the asymptotic 
from these specific points as soon as the distance z from 
the entrance plane exceeds V z cto [or 0, for a > n/2]. 
(Actually, the Laplace image singularities at p = T s and 
p = p n determine the entrance effect, i.e., the waves pro- 
duced by the seed at the entrance aperture amplified in 



the pump field in plasma; these waves will inevitably drop 
behind the laser pulse and their effect therefore vanishes 
at z — > 00.) The contribution from p = 0, 

- < v 1 L (l-r)exp[(V z /V x )r s x} \ 

a s+ (x, z -> 00, s — <^ 1 - V , 

sT s { 1 -rexp[(V z /V x )T s L x \ J 

thus determines the long-scale evolution of the instability 
in plasma, which is dominated by the lateral reflections 
of scattered EMW. Neither s — Cj (where T s — 0) nor 
s = are the singularities of the image a s+ (x, z —* 00, s), 
so contributions to the asymptotic originate from the sin- 
gular points s n only, which are the solutions of the equa- 
tion r s „ = v n [here, v n = {V x /V z ){- lnr + 2nin)/L x ], or 
s n + VnVzS^ — iG 3 = 0, n is integer. The fundamental 
specific point sq with the maximum real part is the root 
of the cubic equation s 3 + vq V z s — iG 3 = which, for 
arbitrary Vq, admits quite a cumbersome explicit expres- 
sion. We address to the physically interesting limit of 
"low leakage", i.e., vo < 3G/V Z , or r > exp(~3GL x /V x ), 
which means that the scattered EMW gains more energy 
between two reflections than loses due to leakage through 
the wall at one reflection. In this case the lateral con- 
vection is mostly suppressed. Solution obtained via the 
perturbation approach reads so ~ (i + V3)G/2 — voV z /3. 
Contribution to the asymptotic from that point is of the 
order of e s °^, which grows in time with an increment 
Re so- However, the absolute value of v n grows with n, 
so an evaluation has to be done of the contribution from 
the points s n with large n. We again use the perturba- 
tion approach with the small parameter ji n — G/(\v n \V z ), 
that is, n > GL x /(2irV x ), and find the solutions s n ~ 
—v n V z - iGfil, Sn' 3) ~ ±G v / -(signn)/i n . Obviously, 
Kesj <C G for fj, n <C 1, so that contribution from these 
points to the asymptotic is negligible compared to that 
from the points s„ with n < GL x /(2itV x ), which is of the 
order of e s °^. Expanding the image a s+ (x, z — > 00, s) in 
the vicinity of s — sq, we arrive at the expression asymp- 
totically valid for G£ 3> 1, 

a * +M ~ I (inv) l ~ir exp ( So€ - lnr £) ' 

which represents the unstable EMW as a quasi-lD expo- 
nentially growing mode. 
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